The stability of thin-shell wormholes with a phantom-like equation of state 
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This paper discusses the stability to linearized radial perturbations of spherically symmetric thin- 
shell wormholes with a "phantom-like" equation of state for the exotic matter at the throat: V — ujo, 
uj < 0, where a is the energy-density of the shell and V the lateral pressure. This equation is 
analogous to the generalized Chaplygin-gas equation of state used by E.F. Eiroa. The analysis, 
which differs from Eiroa's in its basic approach, is carried out for wormholes constructed from the 
following spacetimes: Schwarzschild, de Sitter and anti de Sitter, Reissner-Nordstrom, and regular 
charged black-hole spacetimes, followed by black holes in dilaton and generalized dilaton-axion 
gravity. 
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I. INTRODUCTION 

A powerful theoretical method for describing or math- 
ematically constructing a class of spherically symmetric 
wormholes from black-hole spacetimes was proposed by 
Visser in 1989 1]. This type of wormhole, constructed 
by the so-called cut-and-paste technique, is commonly 
known as a thin-shell wormhole, since the construction 
calls for grafting two black-hole spacetimes together. The 
junction surface is a three-dimensional thin shell. The 
cut-and-paste technique is now considered standard. 

While there had already been a number of forerunners, 
the concept of a traversable wormhole was proposed by 
Morris and Thorne in 1988 0. Ten years later a renewed 
interest was sparked by the discovery that our Universe is 
undergoing an accelerated expansion 0,13]: a > in the 
Friedmann equation a/a = — ^-(p + 3p). (Our units are 
taken to be those in which G — c — 1.) The acceleration 
is caused by a negative pressure dark energy with equa- 
tion of state (EoS) p = ujp, ui < — I, and p > 0. A value 



of uj < — i is required for an accelerated expansion, while 
uj = — 1 corresponds to a cosmological constant ]5]. The 
case uj < — 1 is referred to as phantom energy and leads 
to a violation of the null energy condition, a primary 
prerequisite for the existence of wormholes. Wormholes 
may also be supported by a generalized Chaplygin gas [(| 
whose EoS is p = —A/p a , where A > and < a < 1. 

In a thin-shell wormhole the exotic matter is confined 
to the thin shell. This suggests assigning an equation of 
state to the exotic matter on the shell. Eiroa [7[ used 
the above generalized Chaplygin EoS V — Aj\o\ a , where 
a is the (negative) energy-density of the shell and V the 
lateral presssure, to perform a stability analysis for lin- 
earized radial perturbations. In this paper we will con- 
sider, analogously, the EoS V = ujo-, uj < 0, which will 
be called a phantom-like equation of state. The stabil- 



ity analysis will be carried out for several spacetimes: 
Schwarzschild, de Sitter and anti de Sitter, Reissner- 
Nordstrom, and regular charged black hole spacetimes, 
as well as black holes in dilaton and generalized dilaton- 
axion gravity. The phantom-like equation of state yields 
explicit closed-form expressions for a. Our approach to 
the stability analysis is therefore different from Eiroa's. 



II. THIN-SHELL WORMHOLE 
CONSTRUCTION 

Our starting point is the spherically symmetric metric Q 
ds 2 = -f(r)dt 2 + [f(r)}' 1 dr 2 +h(r)(d9 2 +sm 2 9d(f) 2 ), (1) 



where f(r) and h(r) are positive functions of r and h(r)is 



increasing. (In Sections III- VI, h(r) 



As in Ref. 



the construction begins with two copies of a black-hole 
spacetime and removing from each the four-dimensional 
region 

n± = {r<a\a>r h }, (2) 

where r = r^ is the (outer) event horizon of the black 
hole. Now identify (in the sense of topology) the time- 
like hypersurfaces 



{r = a | a > r h }. 



The resulting manifold is geodesically complete and pos- 
sesses two asymptotically flat regions connected by a 
throat. Next, we use the Lanczos equations [H0-[IH 



S 1 



"ft? ^ 



(3) 



where [Kij] — Kf- — K^- and [K] is the trace of K % j. 
In terms of the surface energy-density a and the surface 
pressure V, 5^ = diag(— a, V, V). The Lanczos equations 
now yield 
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and 



V = - {[K\] + [K%]) 



(5) 



A dynamic analysis can be obtained by letting the ra- 
dius r — a be a function of time Hi. As a result, 



1 

2ira 



and 



1 2a + /'(a) 



(6) 



(7) 



Since a is negative on the shell, we are dealing with exotic 
matter. In fact, the weak energy condition (WEC) is 
obviously violated since the radial pressure p is zero for 
a thin shell. (The WEC requires the stress-energy tensor 
T a p to obey T a p\x a pP > for all time-like vectors and, 
by continuity, all null vectors.) So for the radial outgoing 
null vector (1, 1, 0, 0), we therefore have T a ^ix a pP = p + 
p = a + < 0. 



III. SCHWARZSCHILD WORMHOLES 

For our first case, the Schwarzschild spacetime, h(r) = r 2 
in line element (J]), as noted earlier. Also, recall that the 
radius r — a is a function of time. It is easy to check that 
V and a obey the conservation equation 

-(va 2 ) + V-(a 2 ) = 0. 

(In Eqs. © and ([7]), the overdot denotes the derivative 
with respect to r.) The equation can be written in the 
form 



da 2, 

— + -(<j + V)=0. 
da a 



(8) 



For a static configuration of radius a , we have a = 
and d — 0. Moreover, we will consider linearized fluc- 
tuations around a static solution characterized by the 
constants ao, ao, and Vq. Given the EoS V = uia, Eq. 
([8|) can be solved by separation of variables to yield 



/a \ 2("+i) 

|ffo1 (t) 



where cto = er(ao). So the solution is 

er(a) = cr J , (To = c(ao). (9) 

Next, we rearrange Eq. ([6]) to obtain the equation of 
motion 



Here the potential V(a) is defined as 

V(a) = f(a) - [2naa(a)} 2 . (10) 
Expanding V{a) around ao, we obtain 

V(a) = V(a Q ) + V'(a Q )(a-ao) + ^V"(a )(a-a ) 2 



+O[(a~a ) 3 ] 



(11) 



Since we are linearizing around a ~ ao, we require that 
V(ao) = and V'(ao) = 0. The configuration is in stable 
equilibrium if V"(a ) > 0. 

Now recall that for the Schwarzschild spacetime, 
f(r) = 1 - 2M/r. It follows that 



2M , , , 
l/(a) = 1 47r 2 aV 



-, 2M ,222 f a o\ 

= 1 47r a en — 

a V a / 

from Eq. ©. From Eq. © with a = 0, 
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1 / 2M 
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so that 



V(a) = 1 



2M 
a 



2M \ a 2 +^ 
—J ^2+iJ- 



(12) 



The first requirement, V^(ao) = 0, is clearly met, but not 
the second. (If the exotic matter on the shell were not 
required to meet the extra condition in the form of an 
EoS, then V'(ao) would indeed be zero @.) From 



V'(a ) = 
2M 



2M , , | | , ;, , . 

1 a (-2-4w)a = 0, 

ao 



we obtain the condition 



1 a /M - 1 
2a /M-2' 



(13) 



Observe that as ao — > +oo, u — > —1/2—, and as ao — > 
2M+, ui —> — oo. At a = 3M, a; = —1. Substituting in 



y"(a) 



4M 


/ _2M 






a 3 





ao +4ti, (2 + 4w)(3 + 4w)a" 



-4 — 4lu 



and simplifying, we obtain the intermediate result 



V"[ao) = -5 



1 a /M - 4 



g V a o/M ao/M a /M - 2 



>0. (14) 



Since the Schwarzschild black hole has an event horizon 
at r = 2Af, a^/M — 2 > 0, and we conclude that the 
inequality V"(a ) > can only be satisfied if 



d 2 + V» =0. 



a < 0. 
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As a result, there are no stable solutions. 

To allow a comparison to some of the other cases, let 
us choose (arbitrarily) a /M = 5, as a result of which 
u = —2/3, and plot V(a) against a/M, as shown in Fig. 
1. 
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FIG. 1: The wormhole is unstable. 

The more general analysis in Ref. [8[ depends on 
the parameter (3 2 (a) = 8V /da, where (3 is usually in- 
terpreted as the speed of sound, so that < 1 < 1. 
There are no stable solutions in this range. However, as 
discussed in Ref. 0, since we are dealing with exotic 
matter, this assumption may be questioned, that is, /3 2 
may be just a convenient parameter. In that case, some 
stable configurations may not be out of question. Our 
additional assumption, the EoS V = log on the shell, 
eliminates this possibility. 



IV. WORMHOLES WITH A COSMOLOGICAL 
CONSTANT 



Observe that V(ao) = 0. As before, we have to determine 
the condition on w so that V'(ao) = 0: 



1 1 - l/(a /M) - {2/ 3) AM 2 (ao/M) 2 
2 1 - 2/(a /M) - (l/3)AM 2 (a /A/) : 



(16) 



(As in the Schwarzschild case, as ao — > +oo, to — > —1/2—, 
and ui — !• — oo as ao approaches the outer event horizon.) 
Substituting in V"(ao) and simplifying, we get 



V"(a ) = 

2 -l/(a /M) 



■ 3AM 2 (a /M) - (2/3)Atf 2 (a /M) 2 



1 



2/(a /M) - (l/3)AA/ 2 (a /M) 2 

> 0. 



(17) 



The form of V"(ao) forces us to consider two cases, a 
positive and negative denominator. 
If the denominator is positive, then 



AM > 



1 



{a /M)[3(a /M) - {2/Z)(a /Mf 



(18) 



This inequality implies that ao/M < 4.5 to keep the right 
side positive. It is easy to show analytically that AM 2 > 
1/9; in fact, (3,1/9) is a minimum. We may also plot 
AM 2 against a/M, as shown in Fig. 2. So for this case, 




A. Schwarzschild-de Sitter spacetimes 



FIG. 2: K = AM 2 is plotted against a/M. 



In the presence of a cosmological constant, f(r) = 1 — 
^ - \Ar 2 . For the de Sitter case, A > 0. To keep 
f(r) from becoming negative, we must have AM 2 < 1/9. 
This condition results in two event horizons, where the 
inner horizon is between 2M and 3M. (See Ref. for 
details.) We therefore assume that a is greater than the 
outer horizon. Proceeding as in Sec. IIII1 



V(a) = 1 - — 



2M 
a 



-Aa 2 - 



«o 



lay 



(15) 



the condition V"(ao) > cannot be met (since we must 
have AAf 2 < 1/9), and we get only unstable solutions. 
Plotting V(a) around ao/M = 5 yields a graph that is 
very similar to the graph in Fig. 1. 
For the second case, 



1 - 



a /M 3 \MJ 



< 



in inequality (1171) we obtain 



AM 



3 - 1 (^X 
\MJ 3 \mJ 



< 



a /M ' 



(19) 



(20) 
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If ao/M > 4.5, then the left side is negative, and the 
condition is automatically satisfied. If ao/M < 4.5, then, 
according to Fig. 2, AM 2 < 1/9, the region below the 
graph. So we conclude that in the second case, the worm- 
holes are stable. 

For comparison, let us choose ao/M — 5 again and 
AM 2 = 0.11 < 1/9, resulting in u = -1.63. The plot of 
V(a) against a/M is shown in Fig. 3. 



So the wormhole is stable whenever 



v ' (a /M)[3( a o/M)-(2/3)( ao /M) 2 ] 

and 

(2) a /M > 4.5. 
V. REISSNER-NORDSTROM WORMHOLES 




a/M 



FIG. 3: The wormhole is stable. 

In summary, in the Schwarzschild-dc Sitter case, the 
thin-shell wormholes are stable if, and only if, 



1 



ao/M 3 



iAM*(^<0. 



B. Schwarzschild-anti de Sitter spacetimes 

To study the case A < 0, we return to inequality (|T7|) 
and consider first a negative denominator: 



1 - 



ao/M 3 \M) 



< 0. 



Solving for AM 2 , we obtain 

AM2 > 3-6/( ao /M) 
(ao/M) 2 ■ 

Since ao/M > 2, we conclude that AM 2 > 0, so that this 
case cannot arise. 

Reversing the sense of the inequality, we have from 
inequlity ([T7| 



AM 



of a o\ _ 2 f a o\ 

6 Km) 3\m) 



> 



i 



ao/M' 



Then the second factor on the left must be negative, 
which implies that ao/M > 4.5. 



If the starting point is a Reissner-Nordstrom spacetime, 
then 



/ ( r ) = 1 ~ — + 72"' 



(21) 



where M and Q are the mass and charge, respectively, 
of the black hole. For < \Q\ < M, this black hole has 
two event horizons at r = M ± y/M 2 — Q 2 . As usual, we 
require that r = a is larger than the outer horizon. 
Here we have 

2M Q 2 
V(a) = 1 + \ 



a a 

Once again, V^(ao) = 0. From V'(ao) — we obtain 



(22) 



1 



{ao/M) 2 - a /M 



2 (a /M) 2 - 2(a /M) + Q 2 /M 2 ' 



(23) 



Substituting into V"{ao) and simplifying, yields the 
following inequality: 



V"(a ) = 

2 -a /M - (Q 2 /M 2 )[l/(a /M)} + 2Q 2 /M 2 



a 2 (ao/M) 2 - 2(a Q /M) + Q 2 /M 2 



> 0. 



(24) 



Since ao/M > 2, the denominator is positive. Solving for 
Q 2 /M 2 , leads to 



> 



a /M 



\Q\ 

M " v /2(a /M) - 1 ' 



(25) 



which exceeds unity. To meet this condition, \Q\ would 
have to exceed M. 

So to obtain a stable solution, we will have to tolerate 
a naked singularity at r — 0, but since ao > 0, the naked 
singularity is removed from the wormhole spacetime. 



VI. WORMHOLES FROM REGULAR 
CHARGED BLACK HOLES 

Thin-shell wormholes from regular charged black holes, 
due to Ayon-Beato and Garcia [l6| , are discussed in Ref. 
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[17]. For this black hole 
f(r) = 1 



2M 
r 



2M ( Q 2 

tanh 

r \2Mr 



(26) 



Again, M and Q are the mass and charge, respectively. 
It is shown in Ref. [lj| that the black hole has two event 
horizons whenever \Q\ < 1.05M. Consider next 



V(a) = 1 
1 - 



2M / Q 2 

tanh — — — 

a \2Ma 



2M 
a 

2M 2M , / Q 2 

1 tanh — ; — 

a,Q ao \ 2Mao 



(t)' 



As before, V(a ) = 0, and from V'(oo) = 0, we get 



1 

U= 2 
where 

gM = 

- 1 



-1 



a /M -2 + 2 tanh[Q 2 /(2Ma )] 



(27) 



(28) 



tanh 



( Q 2 



2 /M< 



\2Ma J 2a /M 



sech 



Q 2 



2Ma 



Based on the graphical output, we get only unstable 
solutions. For example, choosing a$/M = 5 again for 
comparison and letting \Q\/M = 0.9, we get uj = —0.63. 
The resulting graph, shown in Fig. 4, resembles Fig. 1. 
Other choices of the parameters lead to similar results. 




a/M 



FIG. 4: The wormhole is unstable. 



VII. WORMHOLES FROM BLACK HOLES IN 
DILATON AND DILATON-AXION GRAVITY 

Of the remaining thin-shell wormholes, based on dila- 
ton and dilaton-axion black holes, respectively, we will 



consider in detail only the latter, which is the more com- 
plicated of the two. 

The dilaton-axion black-hole solution, inspired by low- 
energy string theory, was discovered by Sur, et ai, 18 1 
and is also discussed in Ref. |19| . We need to list certain 
parameters in order to define V(a). As in the Reissner- 
Nordstrom wormhole, there are two event horizons, de- 
noted by r_ and r + , respectively. Returnin g no w to line 
element Q]), we can list both f(r) and h(r) jl8| : 



(r — r_)(r — r + ) 
(r - r ) 2 - 2n (r + r ) 2n ' 

(r + r ) 2n 



h(r) = 



(r - ro) 



2n-2 ' 



Since h{r) is no longer equal to r 2 in line element (fTJ, 
Eq. ([6j) becomes 



1 I ft \ _j_ ' 2 



(29) 



and the conservation equation (jHJ has to be replaced by 
0: 



— (aA)+V — 
dT or 



{[h\a)] 2 -2h{a)h'\a)} ^ f ^ h \ (30) 

where A = Anh(a) is the area of the throat by Eq. ((T|). 
The prime and dot denote, respectively, the derivatives 
with respect to a and r. Substituting Eq. (|2T)1) on the 
right-hand side, we get 

^-[4irh(a)<T}+V^-[4irh(a)} 
dr ' dr 

= -{[ h '(a)] 2 ~2h(a)h"(a)}^l, 



whence 



jL[ah(a)]+V±[h(a)] 
da da 



- { [h'{a)] 2 -2h{a)h"{a)} — 



2h'{a) 



Our final form is 

h{a)a +h'(a)(a + V) 



+ {[h'(a)] 2 -2h(a)h"(a)} — 



2h'(a) 



(31) 



Making use of V — ua, this equation can be solved by 
separation of variables: 



a (a) = <tq 



h(a )' 


3/2+w 


\h'(a o y 


-l 


. H a ) _ 




h'(a) _ 





(32) 
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(Here we used the fact that h'(a) > 0.) It is shown in 
Ref. [HI that 

4[a + (1 - 2n)r ](a - r_)(a - r+) 
°"o = — Z~X7Z — — ~ ) Wo) 



D(a - r )(a + r a ) 



where 



D = 8 7 r(a ^r ) 1 -"(a +r ) V(«o - r_)(a - (34) 

Using the equation of motion d 2 + V(o) = once again, 
we get from Eq. 



V(a) = /(a) - 
Eq. ([32l) now yields 



(35) 



V(a) = 



(a - r-)(a - r + ) 
(a - r ) 2 - 2 "(a + r ) 2n 



A Ha) 



h(a ) 



h'(a) 



3+2u 



(36) 



While it is easy enough to check that V(a ) = 0, it is 
no longer convenient to compute w as a function of the 
various parameters. Plotting V(a) against a instead of 
a/M, we can determine ui by trial and error: V(a) must 
be tangent to the a-axis at a = a , where V(a ) = 
automatically. For example, if ao = 5, ro = 1, r_ = 2, 
r+ = 2.05, and n = 0.8, then w = -0.915. If a = 5, 
ro = 1, r_ = 2, r+ = 3, and n = 0.8, then cj = —1.132. 
Reducing n to 0.6 produces ui = —0.84 in the first case 
and to — —1.041 in the second. In all cases the graphs 
are concave down at ao = 5 and look similar to the graph 
in Fig. 1. Based on the graphical output, there do not 
appear to be any stable solutions. 
For the dilaton case we have [2fJ 



V(a) 



A 



1 - 



B 



(l-6 2 )/(l+6 2 ) 



, Ha) 
h (a) 



h(a ) 
h(a) 



3+2u 



(37) 



where h(a) = a 2 (l — £>/a) 262 /( 1+h2 ' for various constants. 
Once again, one can readily check that V(ao) = 0. 

As in the dilaton-axion case, uj can be found by trial 
and error. For example, if ao = 5, b = 0.5, A = 2, and 
5 = 1, then u = -0.693; if a = 6, b = 0.8, A = 4, and 
B = 2, then ui = —0.94, etc. The resulting graphs are 
similar to those in the dilaton-axion case. 



VIII. CONCLUSION 

This paper discusses the stability to linearized radial per- 
turbations of spherically symmetric thin-shell wormholes 
with the equation of state V = ucr, uj < 0, for the ex- 
otic matter at the throat. This EoS is referred to as 
phantom-like. Various spacetimes were considered. 

It was found that the wormholes are unstable if con- 
structed from Schwarzschild spacetimes, as well as from 
black holes in dilaton and dilaton-axion gravity. For the 
Reissner-Nordstrom case, stable solutions exist only if 



> 



a /M 



\Q\ 

M " y/2(ao/M)-l' 

leading to a naked singularity. For the Schwarzschild-de 
Sitter case, the wormholes are stable if, and only if, 



1 - 



ao/M 3 



3 \MJ 



< 0. 



In the Schwarzschild-anti de Sitter case, the configura- 
tions are stable whenever 



(1) AM 2 < 



1 



(a /M)[3(a /M) - (2/3)(a /M) 2 ] 



and 



(2) a /M > 4.5. 
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